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TAMAGAWA NUMBERS OF DIAGONAL CUBIC SURFACES 
OF HIGHER RANK 

EMMANUEL PEYRE AND YURI TSCHINKEL 
Abstract. We consider diagonal cubic surfaces defined by an equation of the form 



■by" 



-dr 



:0. 



Numerically, one can find all rational points of height SC B for B in the range of up to 
10^, thanks to a program due to D. J. Bernstein. On the other hand, there are precise con- 
jectures concerning the constants in the asymptotics of rational points of bounded height 
due to Manin, Batyrev and the authors. Changing the coefficients one can obtain cubic 
surfaces with rank of the Picard group varying between 1 and 4. We check that numer- 
ical data are compatible with the above conjectures. In a previous paper we considered 
cubic surfaces with Picai'd groups of rank one with or without Brauer-Manin obstruction 
to weak approximation. In this paper, we test the conjectures for diagonal cubic surfaces 
with Picard groups of higher rank. 
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1. Introduction 

This paper is devoted to numerical tests of a refined version of a c onjec t ure of Manin 
ab out th e number of points of bounded height on Fano varieties (see [ BM |, [FMl], [Q, 
or [BT| for a description of the conjectures). The choice of diagonal cubic surfaces to test 
these conjectures was motivated by the work of Heath-Brown [ |H-B | in which he treated 
the cases 



x^ + r-^ 



Z'^ + aT^ 







for a = 2 or 3. The results he obtained were used as a bench mark for the subsequent 
attempts to interpret the asymptotic constants (see, in particular, |S-D], [Pq] and [PT]). 
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2 EMMANUEL PEYRE AND YURI TSCHINKEL 

More precisely, we consider a diagonal cubic surface V C Pq given by an equation of 
the form 

aX^ + bY^ + cZ^ + dT^ = 0. 

Let H be the height function on P''(Q) defined by the formula: for any Q — {x : y : z : t) 
in P'^(Q), one has 

i/(Q) ^max{N,|yUzUi|} if I ^"^'^'"'^^ ^^'' 

[gcd{x, y,z,t) = 1. 

Let U be the complement in V to the 27 lines. We are interested in the asymptotic behavior 
of the cardinal 

NuA^) = tt{ Q e U{Q) I H{Q) ^ B } 

as B goes to infinity. 

Assume that V(Q) is Zariski dense, which by a result of Segre (see [ Man2 , §29, §30]) 
is equivalent to V{Q) ^ 0. It is expected that 

N^^^iB) = BP{log{B)) + o{B) 

as B goes to +oo, where P is a polynomial of degree rkPic(y) — 1, with leading coeffi- 
cient 6jj{V). This constant has a conjectural description. The goal is to compute 6jj{V) 
exphcitly in the examples at hand and to compare it with numerical data. Our previous 



paper [ PT | was devoted to surfaces with Picard groups of rank one with or without Brauer- 
Manin obstruction to weak approximation. In this paper, we consider examples with Picard 
groups of higher rank. Note that in these examples the error term 

(N^^HiB) - eH{V)B{logBY^^''^^''^-^)/B{\ogBY^^'<'''>-^) 

is expected to decrease more slowly. Indeed, if rkPic(y) = 1 this error term is expected 
to decrease as l/B' for some e > 0, whereas for higher ranks it should be comparable 
to l/logS. However we observe a very good accordance, which is even more striking 
if we take into account that a polynomial P of degree rkPic(y) — 1 should appear in 
the asymptotics and use a rather straightforward statistical formula to estimate the leading 
coefficient of this polynomial from numerical data. 

The paper is organized as follows: in section g we define Ofj{V). Section g contains 
the description of the Galois action on the geometric Picard group Pic(y). In section Q 
we compute the Euler product corresponding to good reduction places. In section S we 
explain how to compute the local densities at the places of bad reduction. In section^ we 
determine in each case the value of the geometric constant a{V) defined in §Q. Section ^ 
is devoted to the description of statistical tools we used to analyze the numerical data. In 
section we present the results. 

2. Description of the conjectural constant 

In this section we give a short description of the conjectural asymptotic constant for 
heights defined by an adelic metrization of the anticanonical line bundle (see [Pa] for more 



details and [BT] for a discussion in a more general setting). 



Notations 2.1. For any field E, we denote by E an algebraic closure of E. If X is a 
variety over E, then X{E) denotes the set of rational points of X and X the product 
X XspcclE) Speci?. The cohomological Brauer group Br(X) is defined as the etale co- 
homology group H?^{X, G„J. For any A in Br{X), any extension E' of E and any P in 
V{E'), we denote by A{P) the evaluation of A at P. 
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For a number field F we denote by Val(F) the set of places of F and by ValJF) the 
set of finite places. The absolute discriminant of F is denoted by dp. For any place v of 
F, let F^ be the u-adic completion of F. If v is finite, then ^^ is the ring of u-adic integers 
and F^ the residue field. By global class field theory we have an exact sequence 

(2.1) ^ Br(F) ^ Br(Fj ^^^ Q/Z ^ 0. 

In this paragraph, V^ is a smooth projective geometrically integral variety over a number 
field F satisfying the conditions: 

(i) W{V, ey) = for i = 1 or 2, 
(ii) Pic(y) has no torsion, 

(iii) the anticanonical line bundle LOy^ belongs to the interior of the cone of classes of 
effective divisors k^^^iy) C Pic(F) (g)^ R. 



The adelic space ^(^4^,) of V coincides with the product IltjeVaUF) ^(^i;)- ^Y ['-■^ 
lemma 1], for any class A in Br(y), one has a map p^ defined as the composition 

y^Ap) - e„eVai(F)Br(Fj ^^ Q/Z 

(^i;)i;eVal(F) ^^ (^(-^'-u ))t;eVal(F) ■ 

Then one defines 

V{Apf^= fl kcr(p^) C n^^). 

AeBr(V') 



By the exact sequence (2. 1 ), one has the inclusion 



V{F) c V{ApY 



where ^(F) denotes the topological closure of the set of rational points. Conjecturally 
both sets coincide for cubic surfaces. (See also the text of Swinnerton-Dyerin this volume). 



The Brauer-Manin obstruction to weak approximation, described by Manin in [ Man! ] and 
by CoUiot-Thelene and Sansuc in [ ]CTS| ], is the condition 



UJ 



viA^r ^ y^Ap). 

Let us assume that the height i/ on V is defined by an adelic metric (|| • ||^)i,gvai(F) °ii 
y^ . By definition, this means that we consider LOy^ as a line bundle, that the functions 
II '11^ are u-adically continuous metrics o\\LOy^(F^ which for almost all places v are given 
by a smooth model of V , and that the height of a rational point a; of y is given by the 
formula 

^y^ujy\x)-{Q-\, B(x)^ n \\y\\-^ 

i;GVal(F) 

where LOy^[x) is the fiber of LLiy^ at x. 

If w e Val(F) the Haar measure da;^ on F^ is normalized as follows: 

- Ig dx^ = 1 if z; is finite, 

- da;^ is the usual Lebesgue measure if F^ -^^ R, 

- da;„ = dzdz = 2dxdy if F„ ^^ C. 



EMMANUEL PEYRE AND YURI TSCHINKEL 



The metric || • ||^ defines a measure u;^ ^ on the locally compact space V{F^). In local 
w-adic analytic coordinates x^^ ^ . . .x^^ on V{F^) this measure is given by the formula 



d d 

A •••A 



dx^ dx„ 



1 V ' ' ' 71 



If M is a discrete representation of Gsi\{F / F) over Q, then for any finite place p of F, 
the local term of the corresponding Artin L-function is defined as follows: we choose an 
algebraic closure F of F containing F. We get an exact sequence 

l^I^^D^^ Gal(Fp/Fp) ^ 1 

where D is the decomposition group and / the inertia. We denote by Fr a lifting of the 
Frobenius map to D C Gal(F/i^) (which up to conjugation depends only on p), and put 

1 



Lp(s,M) 



Det(l-(HFp)-^FrjM^p)' 

We fix a finite set S of bad places containing the archimedean ones so that V admits a 
smooth projective model ^ over the ring of ^-integers ff g. For any p in Val(F) — S* we 
consider 

ip(s,Pic(F)) =Lp(s,Pic(F)®zQ). 

The corresponding global L-function is given by the Euler product 

Ls{s,Vic{V))^ n Lp(s,Pic(F)) 

peVal(F)-S 

which converges for Re s > 1 and has a meromorphic continuation to C with a pole of 
order t — rkPic(y) at 1. One introduces local convergence factors Ay given by 

^ ^|L„(l,Pic(F))if«eVal(^)-5, 
1 1 otherwise. 

The Weil conjectures (proved by Deligne) imply that the Tamagawa measure 

veVal{F) 

converges on V{Ap) (see [ pq , proposition 2.2.2]). 

Definition 2.2. The Tamagawa measure on V{Ap) corresponding to the adelic metric 

(II • L).,eVai(F) is defined by 

^H = ^dimy 1™(^ - lYLs{s,PiciV)) n K'^H,.- 
y/dp veVa.i{F) 

From the arithmetic standpoint, it seems more natural to integrate o;^ over the clo- 
sure V{F) C V{Ap) (as in the original approach to the Tamagawa number). However, 
computationally, it is easier to work with V{Ap)^'^. Therefore, following a suggestion of 
Salberger, we define here 

Definition 2.3. 

(2.2) T^{V)=u:^{V{Af^) 
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Let Pic(l/)^ be the dual lattice to Pic(y). We denote by dy the corresponding Lebes- 
gue measure on Pic(T^)^ ®^ R and by 

R|VyeA,ff(F), (.T,2;)>0} 



K^{vY = {xe 


Vic[VY (5 


the dual cone of k^f^iy). 




Definition 2.4. We define 




a{V) 


1 


(^-1)! 



,-('^ 



.v) 



dy 



K«{v)'' 



and 



p{V)^W\k,Vic{V)). 
The theoretical constant attached to V and H is defined as 
(2.3) e^iV) = a{VmV)TH{V). 

In the following sections we compute 9fj{V) for various diagonal cubic surfaces. 

3. The Galois module Pic(F) 

The description of this Galois module is based upon the study of the 27 lines of the 
cubic. We fix notations for these lines extending those given by Colliot-Thelene, Kanevsky 
and Sansuc in [ |CTKS| , p. 9]. 

Notations 3.1. From now on 1/ is a diagonal cubic surface V given by an equation of the 
form 



(3.1) 



aX-^ + bY-^ + cZ^ + dT^ 







where a,b,c and d are strictly positive integers with gcd(a, 6, c, d) = 1. Let 

S — {cxD, 3} U {p I p\abc} 

We fix a cubic root a (resp. a' , a") of b/a (resp. c/a, d/a) (which is assumed to be in 
Q if b/a (resp. c/a, d/a) is a cube in Q) and we put 



(3 



a' 



We also consider 

a" 'ifad 

aa' V be ' 



7 



/3' 



7 



a" 


sId 


a 


y-b 


a 


iab 


a'a" 


~\lcd 



and 



and 



P" 



7 



a 
a 



a 



We denote by 6* a primitive third root of one. The 27 lines of the cubic surface (3.1) are 
given by the following equations, where i belongs to Z/3Z: 




jX+9^aY = 0, 
\Z+9^PT = 0. 

X+e^a'Z = 0, 
Y+e'+^P'T = 0. 

X+e'a"T ^ 0, 
Y+e^+^l3"Z = 0. 




jX+e'aY = 0, 
X+ff'a'Z = 0, 

Y+e'+^frr = o. 

X+e'a"T = 0, 
Y+e'/3"Z = 0. 




ix+e'aY = 0, 
\z+e'+'^f3T = 0. 

X+e'a'Z = 0, 
Y+9'(3'T ^ 0. 

X+e'a"T = 0, 
Y+e'+^f3"Z ^ 0. 
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Let K be the field Q(j, a, a' , a"). It is a Galois extension of Q. In the generic case, K is 
an extension of degree 54 with a Galois group isomorphic to 

(Z/3Z)3 X Z/2Z. 

It is generated by the elements c, t, t' and t" characterized by their action on 6, a, a' and 





e 


a 


a' 


a" 


c 


02 


a 


a' 


a" 


T 


e 


0a 


a' 


a" 


t' 


9 


a 


9a' 


a" 


t" 


9 


a 


a' 


0a" 



Their action on the 27 lines is given as follows: for t we have 

L{i)^L"{i + l) Mil) — >M"{i) Nil) — >N"{{) 

(3.2) \ / \ / and \ / 

M'ii) N'{i) 



for r': 
(3.3) 
forr" 
(3.4) 

for c: 



L'{i + 2) 

Lii)^L"ii) 

\ / 
i'(z) 

Lii)^L'{i) 

\ / 
L"{r) 



M{i) — >M"{i + l) N{i) — >N'{i) 

\ / and \ / 

M'{i + 2) N"{i) 



Mil) — >M'ii) 

\ / 
M"ii) 



N{i) — >N"{i + l) 



and 



N'ii + 2) 
L{0) L'(0)--L"(0) M{0)^M'{0) M"(0) 



i(l) L'(l) i"(l) Af(l) M'(l) M"(l) 

IX XI 

i(2) L'(2) L"(2) M(2) M'(2) M"(2) 



(3.5) 



N{0) N'{0) N"{0) 



N{1) N'{1) N"il) 



N{2) N'{2) N"{2). 



To describe the relations between the classes of these divisors in Pic(y) we consider V 
as the blow-up of a plane P|j in six points P^, P^, P^, P^, P^ and Pg. The 27 Unes may 
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be described as the 6 exceptional divisors E^, . . . ,Eq, the 15 strict transforms L^ of the 
projective lines {P^Pj) for 1 ^ i < j ^ Q and the 6 strict transforms of the conies Q^ 
going through all points except P- . Let A be the preimage of a line of P|^ which does not 
contain any of the points P^ , . . . , Pg. Then 



{[A],[E,],[E,],[E,],[E,],[E,],[E^]) 
is a basis of Pic(y) and we have the following relations in Pic(l/): 

[L^J = [A] - [E^] - [E^] for 1 ^ z < J < 6, 
^^■^^ m^2[A]-J2[E^]. 

In the following, we choose the projection of V to P|^ so that we have the equalities: 



(3.7) 



E^ 


= P(0), 


^2^ 


= ^(1), 


^3^ 


= L{2), 


E,- 


= M{^), 


^5^ 


= M(1), 


E, 


= M(2), 


Qi 


= i'(l), 


Qi 


= L'(2), 


Q3 


= P'(0), 


Q, 


= M"(2), 


Q, 


= M"(0), 


Q, 


= M"{1), 


'1,1 - 


= i"(l), 


L2.,3 


= L"i2), 


-^3,1 ■ 


- i"(0), 


"4,5 ' 


= M'(0), 


A5,6 ^ 


-M'(l), 


LeA 


= M'(2), 


'1,4 - 


-A^'(2), 


Liy- 


= iV'(0), 


L,s- 


= 7V'(1), 


"2,4 ' 


= A^"(0), 


L2., 


= 7V"(1), 


^2fi - 


= ^"(2), 


"3,4 - 


-^(1), 


-^3,5 ^ 


= iV(2), 


L,s 


- A^(0). 



Notations 3.2. We consider the etale algebra E^ over Q defined as Q(7) if ad/hc is not 
a cube in Q and as Q(0) x Q otherwise. Similarly, we define the algebra E^ (resp. E^) 
corresponding to 7' (resp. 7") and we put 

E = E^x E^x E^. 

We also consider the following elements of Pic(y) 

el = [i(0)] + [L(l)] + [i(2)], el = [L'(0)] + [L'(l)] + [L'(2)], 

e^ ^ [i"(0)] + [i"(l)] + [P"2], el = [M"(0)] + [M(l)] + [M'(2)], 

e? = [M'm + [M"(l)] + [M(2)], el = [M(0)] + [M'(l)] + [M"(2)], 

e3 = [iV'(o)] + [iV'(l)] + [N'{2)], el = [7V"(0)] + [7V"(1)] + [N"{2)], 
el = [7V(0)] + [N{1)] + [N{2)] 

and the sets 

'^1 = i^oi ^1, 62), 1^2 ~ 1^0' ^ij 62), fo's = {eo, 61, 62), 
and S' = S'j^US'^U S^. 

Lemma 3.3. The sets S-^, S'2 and S"^ are globally invariant under the action ofGal{K/Cl) 
and the etale algebra corresponding to the set S'^ is isomorphic to E^. 
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Proof. The fact that the sets S-,^, S^ and <S^ are globally invariant follows immediately 
from the descriptions (3.2)-( 3.5). The etale algebra F corresponding to a finite Gal(ii'/Q)- 
set ^ may be defined as the algebra 

(i^[^])Gal(K/Q) 

where Kl,"^] is the algebra K'^ and where G&\{K/Q) acts simultaneously on K and ^. 
In the generic case, let us consider 

cr = tt', a' = t't" and a" == t"t. 

Then a sends 7 on 6*7 and acts trivially on 7', 7" and 0. We may describe similarly the 
actions of a' and a" . The action of Gal(ii'/Q) on (fj^ in the generic case is given by the 
table 





4 


el 


el 


c 


eh 


el 


el 


(7 


el 


el 


eh 


<j' 


eh 


el 


el 


a" 


eh 


el 


el 



This implies that if ad/bc is not a cube in Q, then S'-^ is isomorphic to 

Gal(A7Q)/Gal(i^/Q(7)) 
as a Gal(AYQ)-set. Then the corresponding etale algebra is 

(i^[Gal(i^/Q)/Gal(i^/Q(7))])^"^(^/'2) -^ K^^KK/Qh)) ^ q(^) ^ e,. 

Similarly if ad/bc is a cube in Q, then we may decompose S"^ into two orbits and we see 
that the corresponding etale algebra is Q{9) x Q = E^. The proofs for ^2 ^"d S"^ are 
similar. D 



Lemma 3.4. There exists an exact sequence o/Gal(_R'/Q) modules 

^ q2 ^ q,[S] ^ Pic(F) ®z Q ^ 0. 



Proof. By (3.6) and ( |3.7[ ), we have in Pic(l/) the relations 

el^[E^] + [E^] + [E^], 

el = 6[A] - 2[E^] - 2[E,] - 2\E^\ - ZyE^ - 3[i?5 



3[£;« 



e\^i{k\-2{E^\-2{E^\-2{E^\, 

el = 3[A] - [E,] - [E^] - [E^] - 2[E^] + [E,] - 2[E^l 

e\ = 3[A] - {E^\ - {E^\ - {E^\ - 2{E^\ - 2{E^\ + {E^l 
2 



el 



3[A] - [E,] - [E^] - [E^] + [E,] - 2[E^] - 2[E, 
i[K]-i[E^]-[E^]-[E,]-[E^l 



^6J' 



el^?.[K]-i[E^]-[E^]-[E^]-[E,], 
el = 3[A] - -ilE.,] - [E^] - [E,;\ - [E^] 
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which proves that the natural projection from Q[^] to Pic(l^) (giz Q is surjective. Moreover 
one has the relations 

xes^ xeS2 xes^ 
which gives a homomorphism of Gal(i4r/Q)-modules 

Q2 ^ Q[(f] 
and the exact sequence of the lemma. D 

Notations 3.5. For any prime p and any finite field extension F of Q, we consider the local 
factor (p of the function (^p at p which is defined by 

Cf»= n (i-ttFr)-^ 

{DeVal(F)|ii^} 

Let F be an etale algebra over Q and F = Yliei ^i its decomposition in fields. Put 

Cf(s)-n^F.(s) and Ci.,p(s) = n^f.,p(^)- 

For any prime p, we denote by i^pip) the number of components of F iS)q Qp of degree 
one over Q 

Proposition 3.6. With notation as above, for any prime p not in S, one has 

(i)Lp(.,Pic(F)) = ^^, 
(ii) Tr(Frp|Pic(F)) = i/s(p)-2. 
Proof. By lemma [3.4[ we have 



L^{s,Vic{V)) 



.^^^^ L^{sMm 



Lpi-sMr 



Thus it is enough to prove that if E is an arbitrary etale algebra over Q corresponding to a 
Gal(Q/Q)-set <#" and if p is a prime such that E/Q is not ramified at p, then 

This well-known assertion follows from the fact that the components of E (g) Q are in 
bijection with the orbits of Fr in S', and the degree of each component is the length of the 
corresponding orbit. This proves (i). 
But this also shows that 

Tr(Frp | Q[^]) - ^^(p) 

which implies (ii). D 



Remark 3 J. Thus the factor A' which was defined in proposition 5.1 in |PT| coincides 



with L {1, Pic(y)) at the good places (as suggested by the referee of that paper). 
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4. EULER PRODUCT FOR THE GOOD PLACES 



We need to compute the number of solutions of ( 3. 1 ) modulo p for all primes not in S. 
Proposition 4.1. For any prime p not in S, one has 

P^ P P^ 

where E is the etale algebra defined in ^J. 

Proof. By a result of Weil (see [ yian2| , theorem 23. 1]), 

Wi^p) - 1 + Tr(Frp | Pic(F))p + /. 
Proposition 3.6 implies that 

Tr(Frp | Pic(F)) = y^ip) - 2. D 



Remark 4.2. W e co uld have proved this result directly as in [pi|]. Let N{p) be the number 
of solutions of (3J_) in Fp. We have 

By @ §8.7 theorem 5], one has 

N{p) = p^ +^Xi{a)X2ib)Xz{c)Xi{d)JQ{Xi,X2iXi,Xi). 

where the sum is taken over all quadruples {xi , • ■ • , X4) of nontrivial cubic characters from 
F* to C* such that X1X2X3X4 = 1 ^^d where 

4 

Jo{xi,x2^x3,Xi) = XI n^*(**)' 

the characters being extended by Xi (0) — 0. For p = 2 mod 3 there are no nontrivial char- 
acters and the formula is obvious. Otherwise there are exactly two nontrivial conjugated 
characters x and %. By [PT, proof of prop. 4.1], we have 

•^o(Xl,X2.X3>X4) ^P{V- 1) 

and 

tty(Fp) = l+p(l + Xxi(a)X2WX3(c)X4(rf))+p' 
where the sum is taken over the same quadruples as above. The formula 

^Xi{a)X2i^)Xi{c)Xi{d) = 

/ab\ /ab\ /'ac\ f ac\ f ad\ / ad\ 

A-J +n J +4m) ^<m) ^At^) +Hd 

implies the result. 



Notations 4.3. For any place v of Q, we put 

r Ce,„(s) 



1 otherwise. 



J if u is finite, 
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Remark 4.4. By proposition |3.6|, A^ = Lp(l,Pic(y)) if p G Val(Q) - S. Thus the 
Tamagawa measure a;^ is given by the formula 






By lemmata 3.2 and 3.4 in [ pT] ] and lemma 5.4.6 in [|Pq|, for any p in Val(Q) — S one has 



^njViQp)) = 



fl^(F, 



p' 



p2 



(see also [ pe^ lemma 2.2.1] and [|PT|, remark 5.2]). Therefore, the local factor at a good 
place p is given by 

(l - i)'(l + ^ + ^) ifp = 1 mods and z/^(p) = 9 

(l - I) '(l - ^) (l + I + ^) if P ^ 1 mod 3 and z.^(p) - 6 

(l ~ p) (l - ^)'(l + P + ^) ifP ^ 1 mods and z.^(p) == 3 

(l - I) "' (l - ^)'(l - f + f) if P = 1 mod 3 and z.^(p) = 

(i-|)(i-f)'(i + I + f) ifP^2mod3. 

We get (for the good places) the factors Cg, C^, Cj and Cg where 

-0^ n (-^)(-^r. 

p=2 mod 3. 

p/Sabcd, 
j3=l mod 3, 
i/_b(p)=9. 



c.- n (i-i)(i-^)*(i * ' 



p^Sabcd, 
p=l mod 3, 



^3- n 0-^) 



p-^ / \ p/ \ p p 



1 \3 



p=l mod 3, 

These products converge rapidly and are easily approximated. 

5. Density at the bad places 
In this section we restrict to cubic surfaces with equations of the form 

(5.1) X^ + Y^ + qZ^ + q^T^ = 
with q prime and 

(5.2) aX^ + aY^ + qZ^ + qT^ = 
with q prime and a an integer coprime to q. 



12 
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Notations 5.1. If y is defined by the equation (3.1), and p is a prime, then we consider 

N*{p'') = jl{(a;,y,z,t) G {Z/p''Zf~{pZ/p''Zf \ ax^+by^+czJ+dt^=0 in Z/p''Z} 



Remark 5.2. By [PT, lemmata 3.2 and 3.4], there is an explicit integer r^ such that 



^H^^iQp)) 



N*{p''o) 
p'^^o 



1-p-l 

If p = 3 and 3 / ahcd, then a direct computation in (Z/9Z)^ gives the value of N*{^) 
and thus of lj^j p(V^(Q„))- Thus, in the following lemma we restrict to the case when V is 
given by (5.1) or (5.2) andp = q. 



Lemma 5.3. If V is given by the equation 



X' 



2rTi3 



then for r ^ 2, 



N*{p'- 



Y^ + pZ^ + p^T 



1- 



3r 



F 

If V is given by the equation 

with p I a, then for r ^ S, 

N*{p'' 



ifp = 2 mod 3, 
3 1 1 — - ) ifp = 1 mod 3, 
I ifp = 3. 



aX^ + aY'^ + pZ'^ 



V 



,3r 




ifp = 2 mod 3, 
ifp = 1 mod 3, 
ifp = 3. 



Remark 5.4. This lemma implies that if V is given by the first equation then the local factor 
at p is given by 

[(i-f)(i-|) ifP^2mod3, 
\'^i/,p(^(Qp)) - < 3(l-i)^ ifp =1 mod 3, 

i| ifp = 3, 

and if V is given by the second equation then this factor is 

[(l-^)' ifp =2 mod 3, 

\'^ff,p(^(Qp)) = < sfl- i] (l-^) ifp =1 mod 3, 

ifp = 3. 



16 
27 



Proof Let us consider the set of quadruples {x, y, z, t) in (Z/p''Z)^ — {pTa/p^ZY such 
that 



(5.3) 



x'^ + y'^ + pz'^ + p^r = in Z/p'^Z. 



If p I a; then p\y, p\z and p | i. Therefore, for any {x,y,z,t) as above, p / a; and p K y. 
But for any triple (y, z, i) in (Z/p''Z — pZ/p''Z) x (Z/p^Z)'^, there exists exactly one a; 
verifying (5.3) if p = 2 mod 3 and exactly three of them if p = 1 mod 3. If p = 3 and 
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y belongs to Ti/yX — iTi/^Ti then ( |5.3[ ) implies that 3 1 z. For any triple {y, z, t) with y 
in (Z/3''Z) - (3Z/3''Z), z in (3Z/3'^Z) and t in (Z/S'^Z) there exists exactly three x in 
Z/3''Z which satisfy (O). We get that 



(p— l)p'" xp'"xp'" 

^*(P'') ^ ) q (p-1)p'"^xp''xp'' 
r,3r 1 -^ p3'- 

q 2x3''~^x3''~^x3'' 



pa 



1-i 

p 




ifp = 2 mod 3, 


-0- 


-J) 


if p = 1 mod 3, 


2 
3 




ifp==3. 



3^ 

Let us now turn to the set of (x, y, z, t) in {Ti/p^TiY — (pZ/p'"Z)'* such that 

ax^ + a?/^ + pz^ + pi^ = 0. 
We decompose this set as follows 

Nlip^) = tt I (x,y,z,t) e (Z//Z)4-(pZ/p'-Z)4 
As above we have the formula 

ivr(p'') 



I ax^ + m/ + pz"^ + pt^ = 0. 

\p\x, p K z, 

I aa;^ + at/ + pz^ + pt^ = 0. 



F 



3r 



(p-l)p'' ^Xp''xp'' _ -. _ 1 
p3r p 

= <3x (p-i)p'-;^xp''xp" ^3(1 1 



)2x3''~^x3''~^x3'' 



P 



if p = 2 mod 3, 
if p = 1 mod 3, 
ifp = 3. 



where for p = 3 we use the equality 

S"--^ X 3"^ - tJ{ (z, t) £ {ZjrZf I z^ = f3 j^qJ 3 1^ 
On the other hand, 



r— lr7\4 



^2{f)^v'{ (x,2;,z,t)e(Z/p'-^Z) 



I ap^x^ + ap^y^ + z^ + t^ == 0. 



and 



iV*(p'-) _ p2 



(p-l)p'' ^xp'' ^xp'' _ -, _ 2_ 



po 



^3. -^><P ''"''^.-'-' ^'''' =\^--v) if^-2mod3, 



q 2x3''~^x3''~^x3''~^ _ 9 
^ 33('-i) ~ ^ 



if p = 2 mod 3, 
ifp= 2r 
ifp = 3. 



We conclude: 



T)3r 



l--!- + -^-4 = l-^ ifp = 2mod3, 



2 I 2 _ 4 
,3^3 3 



if p = 1 mod 3, 
ifp = 3. D 



6. The constant a{V) 

Since the cubic surfaces we consider in this paper are Q-rational (which implies that 
P{V) — 1), it remains to compute the rank t of the Picard group and the value of ct{V). 
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Proposition 6.1. IfV is given by the equation 

(6.1) X^ + Y^ + aZ^ + aT'^ = Q, 

where a is not a cube in Q, then rkPic(V^) — 2 and a{V) — 2. 
IfV is given by the equation 

(6.2) aX^ + oY^ + hZ^ + hT^ = 0, 

where a and b are strictly positive integers and b/a is not a cube in Q, then rkPic(F) = 3 
and a{V) — 1. 

IfV is given by the equation 

(6.3) X^ + Y^ + Z^ + T^ = 
then rk Pic(F) = 4 and a{V) = 7/18. 

Proof. To compute a{V) we shall use its original definition [ Pq , §2]: 

aiV) - Vol{x e A^ff(F) I {ujy\x) = 1} 

where the Lebesgue measure on the affine hyperplane 

jr(A) = { X e Pic(F)'' ^z R I {ujy\x) = A } 

is defined by the {t — l)-form dx which is characterized by the relation 

dx A da;^^ = dy 

(where dujy^ is the linear form defined by cjy^ on Pic(F)^ and dy is the form correspond- 
ing to the natural Lebesgue measure on Pic(F)^ ^2, ^)- More explicitely, let (cj^, . . . , ej 
be a basis of Pic(F) and (e^, . . . , e^) be the dual basis. Write 



" = E^^ 



with Aj 7^ 0. Let f^, . . . , f^_^ be the projection of e^, . . . , e^_]^ on J^(0) along e^. Then 

da; =^d/i^ A---Ad/,^_i. 

When V is given by the equation ( |6.1[ ) the Galois group Gal(iir/Q) is 

Z/3Z >^ Z/2Z 

and the orbits of its action on the 27 lines are 

Oi-{L(0),L'(0),L"(0)}, 
02-{i(l),i(2),i'(l),L'(2),L"(l),L"(2)}, 

03 = {M(0),M(2),M'(0),M'(1),M"(1),M"(2)}, 

04 = {A/(l),Af'(2),M"(0)}, 

05 = {NiO),N{l),N'{l),N'{2),N"iO),N"i2)}, 
O,^{N{2),N'{0),N"{l)}. 

In the basis ([A], [£;i], . . . , [^g])' ^basis of Pic(T/) = (Pic V)Gai(i^:Q) is given by 



tj 



= [^4] - 2[i?5] + [^6 



y , 02 — [J^4J - ^[J^5J -r i-L^Qi- 
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The effective cone Agg(F) is generated by the classes [OJ = J2xeo ■^' which in the basis 
(60,6;^) are given by 



[O4] = e, - 



[0,] = 2e„ 
[O5 



2e, 



[O3] = 2e, 



-2' 



Therefore, this cone is generated by the elements e^ — Cj and e^^ + Cj and a{V) is given 
as the volume of the domain 

X = 1, x + y>0 and x — y > 0, 

that is, as the volume of the segment [—1,1] and a{V) — 2. 

If V is given by the equation (5.2) then Ga^X/Q) is isomorphic to 

Z/3Z X Z/2Z 

and the orbits of the Galois action on the 27 lines are 

O, = {L(0)}, 

0, = {L{l),L{2)}, 

03 = {L'(0),L"(0)}, 

04 = {L'(1),L"(2)}, 

05 = {L'(2),L"(1)}, 

Og = {M(0), M(l), M(2), M' (0), M'(l), M'(2)}, 

O7 = {A'f"(0),M"(l),M"(2)}, 

Og = {N{0),N{1),N{2),N"{0),N"{1),N"{2)}, 

O,^{N'{0),N'{l),N'{2)}. 

A basis of Pic(V^) is given by 



m, 



[E2] + [E,], 



and the cone A^ff{V) is generated by 



[Oil =63, 
[O4] = ei + 62 
[O7] = 2ej^ — 65 



-3' 



-3' 

63, 



[O2 

[O5] = ei - 63, 

[Og] = 2ei + 262 - 63 



[O3] = ei - 62, 
[Og] =61 + 62 + 6, 
[O9] = ei - 262 



3' 



that is, by 



'2' 



63: Gj^ + 62 



-3' 



Ze^ 62 63, 



26o 



(since 3[03] = [Oj] + [Og]). Thus a(y) is the volume of the domain given by 

X — 1, y > 0, z > 0, 
X + y — z > 0, 
2x - y - z > 0, 
^x -2y + z > 0. 
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Using the description above, a{V) is the volume of 



< y, < z, 
z~-y <l, 
y + z<2, 
[2y-z<l. 




Therefore a{V) — 1. 

If V is given by the equation (^3]), then Gal{K/Q) = Z/2Z and the orbits of the 
Galois action on the 27 lines are given by 

O, = {L{0)}, O, = {L{1), L{2)}, O, = {L'{Q), L"(0)}, 

0,^{L'{1),L"{2)}, 0,^{L'{2),L"{1)}, 

Oe = {M(0), M'(0)}, O, = {Af (1), M'(2)}, Og = {M(2), M'(l)}, 

Og = {Af"(0)}, Oio = {M"(l),M"(2)}, 

O,, = {7V(0),iV"(0)}, Oi2 = {7V(l),iV"(2)}, Oi3 = {iV(2),iV"(l)}, 

Oi4 = {7V'(0)}, 0,,^{N'{1),N'{2)}. 



A basis of the Picard group is given by 

e, = [A]-[E,], e^ = [E,l e, = [E,] + [E,], e^ = [E^] - 2[E,] 
The effective cone A^g{V) is generated by 



[^6 



[CgJ ~ 6e^ — 2e2 — e^ — e^, 

[O5] = Se^ - 262 ~ 63 - 64, 

[O7] = ei - 64, 

[CgJ = 2ej^ ~ ^2 — 63 — 64, 

i^iiJ ~ ^^1 ~ "^3 ~ ^4' 
[O15] = 2ei - 262 - 64- 



[O2] = 63, 

[O4] = 36^ - 263 - 64, 
[Ow] ^ 46i - 262 - 263 - 64, 



\o 



12J ~ 264 63 64, 



i^l4J ~ ^1 ^2' 



Since [O3] = [O5] = [O9] + [O44] and [O,,] - [O42] = P,] + [O^], we get that A^^{V) 
is generated by 



^21 63, 064 263 



-4, 



64 64, Z64 



-4' 



TzC-i ^^^2 ^oq ^4.1 ^6-1 ^S' 1 2' ^L-i ^Cij 

The anticanonical class is given by 



Wy = 364 



60 60 6^ . 
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Thus a{V) is the volume of the domain 



3x — y — z — t — 1, 
y>0, z>0, 
x-y > 0, 
2x- z > 0, 
x-t>0, 
3x - 2z - t > 0, 
2x — y — z — t>0, 
4x - 2y - 2z - t > 0, 
2x-2y-t >0, 



that is, of the domain P in R'' given by 



y>0, z>0, 
x~y>0, 
2x - z > 0, 
l-2x + y + z>0, 
1 + y-z > 0, 
1-x >0, 
1 + x — y — z>0, 
1 — X — y + z>0. 



We compute its volume as follows: decompose P into cones with appex (0, 0, 0) and 
supported by the faces not containing this point. Thus we consider the following faces of 
P: 



F, 



1 - x = 0, 
l + y-z^Q, 
1— X — y + z — Q. 



F^ : 1 - 2x + y + z = 0, 
F^ : 1 + x — y — z = 0, 



One has 



1 ^ 
,(y)=Vol(P) = -^Area(Fj. 



The area of F-^ is the volume of the domain 



fy>0, z>0, 
l~y>0, 
2- z>0, 
-1 + y + z > 0, 
1 + y-z > 0, 
2-y- z>0, 

{z-y >Q, 
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and we get Area(Fj^) = ^. For F^ we have the equations 

'y>o, 

-l + 2x-y > 0, 
x~y >0, 

i + y>o, 

2 - 2x + 2y > 0, 
l-x >0, 
2-a; >0, 
x-2y> 0. 

We get Area(F2) = ^. For F^ we have the same equations and the same area. For Fj we 
have the equations 

'y>o, 

l + x-y >0, 
x-y>Q, 
-l + x + y>0, 
2-a; >0, 
-x + 2y> 0, 
1 -a; >0, 
2-2y > 0. 

We find Aiea{F^) = 1/8 + 1/24 = 1/6. The face Fg is given by the same equations and 
Area(F5) = 1/6. Finally 



y j 


^ 










/ 






X 












X 



aiV) 



1 /I 4 



3V2 6 
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7. Some statistical formulae 



With the program of D. J. Bernstein [ Be | , given a number B, we can compute the value 
of Njjjj{2'') for 100 < 2'' < B. Thus we get a family of pairs {B^,Njjj^iB^))j^^^^j^, 
where B^ are consecutive powers of 2 for i < N and Bj^ = i?. For any i between 1 and 
iV, let 

x^^\og{B^) and y, = JV^^(i3J/i3,. 

We expect an asymptotic of the form 

N^,,{B) = BP{log{B))+o{B), 

where P is a polynomial of degree t—1 with a dominant coefficient equal to 0^ (V^). Thus 
we look for a polynomial Q of degree i — 1 such that 

N 



L 

is minimal and compare the leading coefficient of Q with 0^(y). 



i=l 
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Notations 7.1. Let R{X,Y) be a polynomial in Q[X, F] and denote by {R{X,Y)) the 
mean value of iRix^,y^))i^^^j^, that is, 



1=1 
If t = 2 the leading coefficient of Q (if it is uniquely defined) is given by 

^ (XY)-(Y)(X) 

If t = 3 the leading coefficient is 

(YX^) {Y){X^) ((^-)-(^)g)K(^x)-(>-)(^)) 

"^2"" 777^. 777^2 ({x^)-{x){x^)Y ■ 

(A ; - (A ; (x2>-(x>2 — 

If i = 4, the leading coefficient is 

/'KY3\ /v\/y3\ ((^")-(-Y)(jf;;'))((yx)-(y)(j^)) I3S 



3 /y6\_ /y3\2_MlhiP£K^£!>li_£ 

7 



/y6\ /y3\2 ji.^ ;-\'^n^", 



with 






(X2) - (X> 



2 



5 = (yx2) - (y)(x2) - i^|_i^K^((yx) - (y)(x)). 

In the next section, we denote by 05?'(^) the leading coefficient A^_^. 
8. Presentation of the results 



We consider only cubic surfaces of the form (6.1), (5.2), or (5.3). By [CTKS, Lemme 
1], the corresponding surface V is Q-rational and, in particular, Br(V") = 0. Thus the 
Brauer-Manin obstruction to weak approximation is void and 

•ueVal(Q) 

Moreover, 

/3(T/) = tti/i(Q,Pic(y)) = l. 



By (2.2) and (2.3), the constant d^{y) may be written as 

Using remark 4.4 we get 

e^{Y) = «(!/) lim(s - l)*+2c^(s) X a;^,^(y(R)) 



3 



n Vi^,p(nQ,))xnQ 

p|3abc(i i— 
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where E is the etale algebra defined in 3.2. The residue of the zeta function could have 



been computed directly (see, for example, [ Zo, chapter 4]), but instead we used PARI. The 



volume at the real place is given by the formula 






where o;^ is the Leray form 



-\x,y,z,t) = — — 3 3 ^-:^dxdydz. 



Decomposing the domain of integration (and using the various expressions of the Leray 
form) it is possible to remove the singularities of this integral which is then easily estimated 
on a computer The factors corresponding to the bad places have been described in section 
raand the constants Cq, C^, C^, and Cg may be computed directly as in section H. 

We considered the following examples: for the cubic surfaces with a Picard group of 
rank 2 we used 



(^l) 


x^ + y3 ^ 2z3 + 4r'^ 


= 0, 


(^2) 


x^ + r3 + 5^3 + 25r3 


= 0, 


(^3) 


^3 ^ y3 ^ 3^3 _^ gy3 


= 0. 


For the rank 3 case: 






(^4) 


^3 ^ y3 ^ 2z^ + 2r^ 


= 0, 


(^5) 


x^ + y-^ + 5Z^ + 5^3 


= 0, 


(^e) 


x^ + y^ + 7Z^ + TT-'' 


= 0, 


(^7) 


2^3 + 2y^ + 3Z^ + iT^ 


= 0, 


and for rank 4: 






(^s) 


X^ + Y^ + Z^ + T^ 


= 0. 
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We draw below the corresponding experimental curves in which we compare the value of 

7V^^^(i?)/(B(logi?)*-i) with 0^(y). 



OniV) 



+ 



+. 



"+++ 



-H-H--^ Oh{v) 



+ 



J I I I I I I I I L 



12 3 4 5 6 7 
^1 



9 10 

log(i?) 



+ 

+ + 



+++-Hh+ 



J I I I I I I I I L 



12 3 4 5 6 7 



9 10 

log(B) 
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OHiV) 



+ 



+ 



+++++++ 



OniV) 



J I I I I I I I I L 



12 3 4 5 6 7 



9 10 

log(B) 







J I I I I I I I I I L 



12 3 4 5 6 7 



9 1011 

log(S) 



OhIV) 



+ 



++ 



"++_ 



"++++-^ 



OniV) 



J I I I I I I I I I I 

0123456789 1011 

log(i?) 



+ 



+ 



+++. 



"+ 



++++_ 



-++-^ 



J I I I I I I I I I L 



0123456789 1011 

log(S) 
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OhIV) 



+ 



+ 



++ 



+ 



+++++++_ 



--^ 



OniV) 



J I I I I I I I I I I 

0123456789 1011 

log(i?) 

^7 



+ 



+ 



+ 

+_ 



"++- 



"+'N-+ 



+ 



J I I I I I I I I I L 



0123456789 1011 

log(B) 
5« 
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We finish with tables of numerical results. The value of 0^h^{V) is obtained from the 
pairs (_Bj, Njj ^{B-)) as described in section]^. We denote by C^;. (1) the limit 

C|;,(l) = lim(s-l)*<|;^(s), 

where i, is the number of components of E^. Note that for the examples with a Picard 
group of rank 2, Cj is equal to 1. 



Surface 


s. 


^2 


^3 


B 


20000 


20000 


20000 


Nu,h{B) 


75984 


49608 


78980 


a{V) 


2 


2 


2 


ad/cb 


2 


5 


3 


Cl,(i) 


8.146241 X IQ-i 


1.163730 


1.017615 


ab/cd 


1/8 


1/125 


1/27 


Cl;.(l) 


6.045998 X IQ-i 


6.045998 X 10-^ 


6.045998 X 10-^ 


ac/bd 


1/2 


1/5 


1/3 


^3(1) 


8.146241 X 10^1 


1.163730 


1.017615 


Xs^niViQ,)) 


4/9 


4/9 


4/9 


Po 


2 


5 




A,,u;^(^(Q,„)) 


3/8 


96/125 




Co 


8.306815 X 10-1 


3.493824 x 10-^ 


3.066383 X 10-^ 


c. 


9.540383 X 10"! 


8.704106 X 10-1 


9.762028 X IQ-i 


C3 


9.893865 x 10"! 


9.906098 X 10-1 


9.892790 X 10-^ 


u;^(F(R)) 


3.255161 


1.360417 


2.221359 


0if(^) 


3.413500 X 10-1 


2.290769 X 10-^ 


3.660885 X 10-^ 


NuMB)/OH(.V)B\ogiB) 


1.123839 


1.093332 


1.089213 


oTiv)/eH{v) 


1.008178 


0.958517 


1.021637 
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For the examples with a Picard group of rank 3, Cg is equal to 1. 



Surface 


s. 


^5 


^6 


S, 


B 


100000 


100000 


100000 


100000 


NuMiB) 


3051198 


1976482 


3420784 


1966160 


a{V) 


1 


1 


1 


1 


ad /be 


1 


1 


1 


1 


a,(i) 


6.045998 X IQ-i 


6.045998 X lO"! 


6.045998 X lO-i 


6.045998 X lO-i 


ab/cd 


1/4 


1/25 


1/49 


4/9 


a,(i) 


8.146241 X IQ-i 


1.163730 


1.265025 


1.028996 


ac/bd 


1 


1 


1 


1 


^3(1) 


6.045998 X IQ-^ 


6.045998 X lO^^ 


6.045998 X lO-i 


6.045998 X lO-i 


A3a;^(y(Q3)) 


16/27 


16/27 


16/27 


16/27 


Pa 


2 


5 


7 


2 


K-^HiV{QJ) 


27/64 


13824/15625 


186624/117649 


27/64 


Co 


8.306815 X IQ-i 


3.493824 x lO^^ 


3.066383 X lO-i 


8.306815 X 10-1 


c\ 


9.540383 X 10"! 


8.704106 X 10-1 


9.297617 X 10-1 


8.196347 X 10-1 


C2 


7.827314 X 10-1 


8.112747 X 10-1 


9.228033 X lO-i 


8.294515 X 10-1 


a;^(t/(R)) 


4.105301 


2.347970 


1.910125 


2.430506 


OniV) 


1.895795 X 10-1 


1.291945 X 10-1 


2.184437 X 10-1 


1.290720 X 10-1 


NuAB)/eH{V)B\og{Bf 


1.214249 


1.154191 


1.181448 


1.149252 


oTm/e^iv) 


0.980969 


1.056178 


1.044971 


0.981831 
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For the last example we have C2 = Cg = 1 and E^= £2 = E^ and we get 



Surface 


'% 


B 


100000 


NuAB) 


12137664 


a{V) 


7/18 


C4(i) 


6.045998 X lO^i 


A3a;^(t/(Q3)) 


16/27 


Co 


3.066383 X IQ-^ 


Ci 


5.129319 X 10^1 


a;^(T/(R)) 


6.121864 


OHiV) 


4.904057 X 10-2 


Nu,H{B)/e^{V)B\og{Bf 


1.621894 


oTm/e^iv) 


1.024630 
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